Moments of the longitudinal proton structure function Fl 
from global data in the Q 2 range 0.75 - 45.0 (GeV/c) 2 
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We present an extraction of the lowest three moments of the proton longitudinal structure function 
Fl from world data between Q 2 = 0.75 and 45 (GeV/c) 2 . The availability of new Fl data at low 
Bjorken x from HERA and at large x from Jefferson Lab allows the first determination of these 
moments over a large Q 2 range, relatively free from uncertainties associated with extrapolations 
into unmeasured regions. The moments are found to be underestimated by leading twist structure 
function parameterizations, especially for the higher moments, suggesting either the presence of 
significant higher twist effects in Fl and/or a larger gluon distribution at high x. 

PACS numbers: 13.60.Hb, 14.20.Dh, 12.38.Qk 



Introduction. The suppression of the longitudinal 
deep-inelastic lepton-proton scattering cross section rel- 
ative to the transverse cross section was an important 
early verification of the quarks' spin-1/2 nature. In fact, 
for a point-like quark the longitudinal structure function 
F L is identically zero. For a composite particle such as 
a proton, Fl is small but finite, and its exact value and 
momentum dependence reflect the quantum interaction 
effects between the proton's quark and gluon (or parton) 
constituents. 

In QCD, one of the novel features of the proton lon- 
gitudinal structure function is its strong sensitivity to 
the nonperturbative initial state distribution of gluons. 
The moments of Fl in particular are related to matrix 
elements of local twist-two operators, which can be com- 
puted directly in lattice QCD. Traditionally, the gluon 
distribution g(x) has been largely determined by study- 
ing the Q 2 evolution of the F^ structure function, which 
at high photon virtualities is dominated by the transverse 
cross section. In recent global fits [H-Q g(x) is further 
constrained at low parton momentum fraction x by jet 
production data in hadronic collisions. 

At large values of x, where the cross sections are 
small, the extraction of the gluon density becomes in- 
creasingly difficult, leading to large uncertainties in g{x) 
at x > 0.3. As a result, the higher moments of Fl, which 
are weighted towards higher values of x, are particularly 
challenging. 

Data on Fl are generally difficult to extract from cross 
section measurements, requiring detailed longitudinal- 
transverse (L/T) separations in which experiments are 
performed at the same x and Q 2 but at different en- 
ergy. Historically, the kinematic range spanned by Fl 
data was therefore rather limited, typically concentrated 
in the small- and intermediate-x regions, whereas a pre- 
cise moment analysis necessitates a broad coverage in x 
at fixed Q 2 . Previous moment analyses consequently re- 
quired recourse to model-dependent estimates of the lon- 
gitudinal to transverse structure function ratio |f|, ren- 



dering a precise evaluation of Fl moments and their un- 
certainties problematic. 

Recently, new data on the proton longitudinal struc- 
ture function have been taken at low x from the HI exper- 
iment Q at HERA, and at large x from Hall C at Jeffer- 
son Lab Q . The latter in particular cover a significant x 
range from Q 2 = 4 (GeV/c) 2 down to Q 2 < 1 (GeV/c) 2 . 
Combined with the previous Fl results, the new data 
allow for the first time a direct extraction of the Q 2 de- 
pendence of several low Fl moments over a large range 
of Q 2 . 

In this Letter, we report the results of such an ex- 
traction, with an accurate determination of the low- 
est three moments of Fl, together with their uncer- 
tainties. Comparison of these to moments computed 
from parametrizations of leading twist parton distribu- 
tion functions (PDFs) [lJ-Q 

then allows one to draw con- 
clusions about the (poorly constrained) gluon distribu- 
tion at large x, or the role of higher twist effects in the 
longitudinal cross section. 

Data Analysis. Our analysis is performed in terms of 
the longitudinal Nachtmann moments, defined as Q 
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where the Nachtmann scaling variable £ = 2x/(l + p), 
with p = y/l + 4M 2 x 2 /Q 2 , and M is the proton mass 
[ll [l0| • The Nachtmann moments are constructed to re- 
move from the data the kinematic dependence on the tar- 
get mass llMl3| , thus allowing a direct comparison with 
the Cornwall-Norton moments calculated from leading 
twist PDFs. At very low Q 2 this treatment may gener- 
ate some residual uncertainty if the contributions from 
the threshold region at £ — s- 1 arc large [lij . 

In this analysis only Fl values extracted from dedi- 
cated, experimental L/T separations of the proton cross 
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FIG. 1: (Color online) The Q 2 and x distribution of Fl data 
sets used in this analysis. The new HI data 0] at very small x 
appear clustered around the vertical axis because of the linear 
x scale. 

section data are used. This constraint is a critical re- 
quirement to avoid the introduction of model dependence 
into the M L n ^ extraction, and to accurately estimate the 
uncertainties on the moments. The utilized proton Fl 
data come from a range of experiments at CERN (EMC 
HU, BCDMS Hf|, NMC fp), SLAC (E140X 0, SLAC 
global [HI), DESY (HI Q), and Jefferson Lab (E94110 
|g], E99118 [20j]). The regions of the Q 2 -x space covered 
by the data sets are shown in Fig. [TJ 

Since much of the data were not taken at fixed Q 2 val- 
ues, the Q 2 bins were chosen to ensure the broad coverage 
in x necessary for a moment extraction. For instance, a 
typical bin at Q 2 = 6.5 (GeV/c) 2 included all data in 
the range of 6 < Q 2 < 7 (GeV/c) 2 . To account for any 
Q 2 dependence in the bin, the data were bin-centered to 
the central Q 2 utilizing a combination of global data fits 
which gave good descriptions of the data over the relevant 
kinematic range. A sample of these data bin-centered in 
Q 2 is shown for several central Q 2 values in Fig. [2] For 
the x integration, the data were then binned in x from 
0.01 to 1, utilizing bins of width Ax = 0.01. 

The data shown in Fig. [2] provide the most comprehen- 
sive kinematical coverage of Fl to date. However, some 
regions of x with sparse data remain, especially at larger 
Q 2 . These gaps were filled by utilizing phenomenological 
fits to calculate the structure function at the center of 
any empty x bin. For data with W 2 > 3.9 (GeV/c) 2 , a 



model obtained by a fit [2l| to world data was used, while 
for W 2 < 3.9 (GcV/c) 2 a fit to the resonance region data 
was employed [22[ . (These models were also used for the 
bin-centering in Q 2 , discussed above.) 

The fit values of Fl were also renormalized bin-by-bin. 
For x > 0.4, where experimental data are abundant, the 
fit value was renormalized by the ratio of the fit to the 
error-weighted mean of the real data points at the start 
and end of each empty interval in x. For intervals with 
x < 0.4, where the x gaps are large and data scarce, the 
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FIG. 2: (Color online) Example plots of the Fl data used 
in this analysis for several Q 2 bins in the range 1.75 — 
10.0 (GeV/c) 2 . 

model was renormalized using the error- weighted mean of 
all real data points up to x = 0.4, to prevent a single data 
point at low x from biasing the renormalization factor. 

Since the longitudinal Nachtmann moment defined in 
Eq. includes contributions from both Fl and F2, the 
entire analysis was repeated for the F2 data of the same 
experiments. After filling the gaps in x for each Q 2 
bin, each structure function was integrated as in Eq. JT]) 
and the separate contributions were summed together to 
generate the longitudinal Nachtmann moments M L n ^ for 
n = 2, 4 and 6. The limits of integration were taken from 
x = 0.01 up to the inelastic or pion production threshold, 
x v = [1 + (m 2 +2¥m T )/Q 2 ] -1 , with the pion mass. 
Above this point, which is below x = 1, only elastic scat- 
tering is possible in the laboratory. Results with and 
without the elastic contribution included are presented 
below. 

A Monte-Carlo procedure was utilized to evaluate the 
errors on the moments. For this purpose, 1000 pseudo- 
data sets were generated for each Q 2 bin and for each 
structure function, by sampling for each data point a 
Gaussian distribution with mean value and width equal 
to the value and total error of the data point. For each 
pseudo-data set, the x coverage gaps were filled using 
the method described above, and the data integrated to 
obtain pseudo-moments. This resulted in two distribu- 
tions, one each for the Fl and F2 contributions to the 
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moment. The value of the Nachtmann moment was de- 
fined as the sum of the mean value of the pseudo-moment 
distributions for the Fl and F 2 components. The statis- 
tical error was defined as the sum in quadrature of the 
standard deviations of each of the two pseudo-moment 
distributions. 

In order to estimate a model dependent (systematic) 
error, the process of filling in the gaps in the data was re- 
peated for three other combinations of models. The first 
two used the ALLM parameterization for F 2 [23| and 
the R1990 model for the longitudinal to transverse ratio 
[24] in the region with W 2 > 3.9 (GeV/ c) 2 , in combina- 
tion with two resonance region fits [22J, [25[ . The third 
combination used the fit to world data from Ref. 
for W 2 > 3.9 (GeV/c) 2 and the resonance region 
for lower W. After generating distributions of pseudo- 
moments for all model combinations, the systematic er- 
ror was defined as the maximum difference between the 
original model combination used to fill in the gaps and 
any of the other three combinations. The systematic er- 
ror only has a significant contribution to the total error 
for the first two Q 2 bins of the n — 2 moment; otherwise, 
the statistical error dominates. Finally, the total error on 
each data point was calculated as the sum in quadrature 
of the statistical and systematic errors. 

Results. The extracted n — 2, 4 and 6 longitudinal 
Nachtmann moments are given in TablelH for Q 2 between 
0.75 and 45 (GeV/c) 2 . The values include the measured 
inelastic contributions as well as the elastic component, 
computed from the global proton form factor parameteri- 
zation in Ref. (2(| . The elastic contribution is significant 
only for the lowest Q 2 bins, and decreases rapidly for 
Q 2 > 2 (GeV/c) 2 . The errors are largely driven by the 
uncertainty on the Fl data. 

The experimental Nachtmann moments are shown in 
Fig. 02 and compared with calculations of the Cornwall- 
Norton moments of F L using the MSTW08 0, ABKM09 
0] and CTEQ- Jefferson Lab (CJ) global PDF param- 
eterizations. 

The MSTW08 fit included data on the F 2 and F L 
structure functions in fixed target experiments and 
HERA collider data on reduced DIS cross sections sat- 
isfying Q 2 > 2 (GeV/c) 2 and W 2 > 15 (GeV/c) 2 . The 
kinematic cuts were imposed to avoid the region where 
higher twist (HT) effects may be significant, thereby ex- 
cluding data at high x. Nuclear corrections in deuterium 
DIS were not included in the fit, although these strongly 
affect the d quark and gluon PDFs, even within the cuts 



used [27] . Jet data in pp collisions were also included, 
constraining the gluon PDF at x < 0.3. 

For the C J analysis @ , a similar data set to that used 
by MSTW08 was fitted, although F L data were not di- 
rectly included. The primary constraint here on the Fl 
structure function, and the large- a; gluon PDF, was there- 
fore from scaling violations in the F 2 data. Significantly, 
the kinematical cuts were relaxed to Q 2 > 1.69 (GeV/c) 2 
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TABLE I: The experimental longitudinal Nachtmann mo- 
ments M L n ^ (scaled by 10~ 3 ) extracted from the data along 
with their statistical errors. The inelastic results for the 
n = 2, 4 and 6 moments are given for each Q 2 bin, with the 
elastic contribution shown only for the four lowest Q 2 bins. 



and W 2 > 3 (GeV/c) 2 , increasing considerably the large- 
x coverage afforded by the high-precision SLAC and Jef- 
ferson Lab data. The less restrictive cuts necessitated in- 
clusion of target mass (TM) and HT contributions, and 
nuclear corrections were incorporated using several mod- 
els of the deuteron. Since the DIS data were limited to 
the F2 structure function, only the leading twist contri- 

(n) 

bution to M L could be computed from these PDFs. 

The ABKM09 fit @ used similar cuts to the CJ anal- 
ysis, and also included nuclear, TM and HT corrections, 
but did not utilize Jefferson Lab data. Fits were per- 
formed to DIS cross sections directly, rather than to the 
extracted F 2 , although the exclusion of jet data from the 
analysis weakens the constraints on the gluon PDF at 
x < 0.3. Furthermore, HT terms were included for both 
F 2 and Fl, allowing calculation of moments up to twist 
4. 

Comparison of the measured moments with the PDF- 
based NLO calculations in Fig. [3] (left) shows a turnover 
of the inelastic moments at Q 2 « 2 — 3 (GeV/c) 2 . This is 
due to the effect of the pion threshold x-^, which decreases 
the limit of integration of the real data at low Q 2 values. 
At about the same Q 2 value, the elastic contribution also 
becomes non-negligible. Theoretical calculations of the 
DIS cross sections which neglect the pion threshold and 
integrate up to x = 1 can therefore be meaningfully com- 
pared to data only for Q 2 > 3 (GeV/c) 2 . 

The leading twist calculations are in generally good 
agreement with the data for Q 2 > 10 (GeV/c) 2 . At lower 
Q 2 the fits underestimate the data, particularly for the 
higher moments where large x plays an increasing im- 
portant role. The disagreement with the low-Q 2 , large- 
n data may reflect the poorly constrained gluon PDF, 
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or possibly large effects from higher order perturbativc 
QCD or or higher twist corrections at large values of x. 

The effect of including higher order terms is illustrated 
in Fig. [3] (right), where the total Nachtmann moments are 
compared with moments calculated from the MSTW08 
PDFs [lj at leading order (LO), next-to- leading order 
(NLO) and next-to- next-to-leading order (NNLO). While 
the LO results generally underestimate the data at low 
Q 2 , the agreement progressively improves with increasing 
order. It is only in the highest (n = 6) moment that any 
discrepancy appears, possibly indicating some underes- 
timation of g(x) at high x. There is a well known and 
large uncertainty on g{x) which can also be observed in 
the substantial differences in the NLO calculations from 
different PDFs, as shown in Fig. (left). 

The role of HT contributions in explaining the missing 
strength at small Q 2 can be explored by comparing the 
ABKM09 fit Q with and without higher twist contribu- 
tions. Inclusion of HT corrections improves the agree- 
ment with data, but overestimates the strength some- 
what, even within the relatively large uncertainty of 
the HT contributions. This remains true even in the 
ABKM09 NNLO fit (not shown), where the LT contri- 
bution increases, albeit more slowly than in the MSTW08 
fit, and the HT terms decrease, leaving the total curve 
stable. Since the ABKM09 fit does not include the recent 
Jefferson Lab data 2il IH, it is not directly constrained 
at lower Q 2 and larger x. Minimizing extrapolation un- 
certainties and precisely studying the interplay of leading 
and higher twist contributions in this region will require 
use of the new data in global fits. 

Conclusions. In summary, we have extracted the low- 
est three Nachtmann moments of the proton longitu- 
dinal structure function over the range Q 2 = 0.75 — 
45.0 (GeV/c) 2 from world Fl data augmented by recent 
small- a; measurements by the HI collaboration at HERA 
and large- a; measurements from Jefferson Lab. Although 
the data coverage in x is not absolute, four different com- 
binations of models have been used to fill gaps in the 
data, allowing the moments to be calculated with a rig- 
orous uncertainty analysis. A reasonable estimate of the 
errors has been obtained by a statistical analysis of the 
moments calculated from a Monte-Carlo procedure for 
sampling the data. 

Comparison of the experimental moments with per- 
turbativc QCD calculations using recent global PDF fits 
reveals a need for cither increased high- a; gluon distri- 
bution, which is possible given its large uncertainty, or 
possibly non-negligible higher twist effects. Discrepan- 
cies between data and the PDFs increase as the moment 
order increases, but are reduced by increasing the order 
of perturbation theory. Neglecting higher twist terms 
or higher order calculations in the global fits may over- 
estimate the extracted gluon distributions at large 
shown by the comparison of ABKM09 and MSTW08 cal- 
culations. 
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FIG. 3: (Color online) Longitudinal Nachtmann moments 
M L for n = 2, 4 and 6, as a function of Q 2 , for the mea- 
sured inelastic (circles) and inelastic plus elastic (squares) 
contributions. The left panels compare the data with calcu- 
lations based on global NLO fits to PDFs from MSTW08 Q, 
ABKM09 [3] (including also HT corrections) and CJ @], while 
the right panels show the comparison with the MSTW08 fits 
for different orders in a s (LO, NLO and NNLO). 

Relatively good constraints on the large- a; gluon can 
already be obtained from the scaling violations of i*2, 
provided that a weak cut on W is considered along with 
TM and HT included, as shown in the CJ fits Q. How- 
ever, inclusion of Fl data in global fits is required to 
provide direct constraints on the gluons, with the recent 
Jefferson Lab data at low Q 2 < 4 GeV 2 , in particular, 
facilitating a considerable extension of reach in x. 
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